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Interactions never lower the ground state kinetic energy of a quantum system. However, at 
nonzero temperature, where the system occupies a thermal distribution of states, interactions can 
reduce the kinetic energy below the noninteracting value. This can be demonstrated from a first 
order weak coupling expansion. Simulations (both variational and restricted path integral Monte 
Carlo) of the electron gas model and dense hydrogen confirm this and show that in contrast to the 
ground state case, at nonzero temperature the population of low momentum states can be increased 
relative to the free Fermi distribution. This effect is not seen in simulations of liquid '^He. 



INTRODUCTION 

It is a common assumption that the addition of inter- 
actions to a noninteracting quantum system will broaden 
the momentum distribution and increase the kinetic en- 
ergy (proportional to the second moment of the momen- 
tum distribution). An intuitive understanding of this 
follows from perhaps the first example encountered in 
studying quantum mechanics, the particle in a box, where 
the kinetic energy scales as the box size . More gen- 
erally the kinetic energy, given by the average curvature 
of the wavefunction, is expected to scale as the "localiza- 
tion length""^. Since interactions lead to a nonuniform 
density, i.e. increased "localization" , the kinetic energy 
is expected to increase. 

For condensed matter scientists, two well ingrained 
many body examples of this, illustrated in Fig. |l|, are 
the broadening of the free Fermi ground state momentum 
distribution due to electron repulsion and the depletion 
of the zero momentum condensate in a strongly interact- 
ing Bose system such as ''He 0]. For the homogeneous 
electron gas, Fig. shows the promotion of low momen- 
tum states to higher momentum and the reduction in 
the step discontinuity at the Fermi surface ||^. Similarly 
Fig. 0b shows that the 100% condensation into the zero 
momentum ground state of a noninteracting Bose system 
is reduced to roughly 8% in ''He with the rest going into 
an almost Gaussian distribution Q . 

The proof of this assumption for systems at zero tem- 
perature is an immediate consequence of the ground state 
variational principle applied to the free particle Hamil- 
tonian, Hq, which states that Kq = (^'ol^/^ol^o) is mini- 
mized by 5*0, the true ground state wavefunction of Hq- 
Using any other wavefunction, such as the ground state, 
^'g, for an interacting system, H = Hq + U, leads to 

K^{-^g\Ho\^g) >Ko, (1) 

demonstrating that the kinetic energy, K, of the inter- 
acting system is never lower than the free particle kinetic 
energy, Kq. 
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FIG. 1: Top figure: Broadening of the ground state elec- 
tron gas momentum distribution at Ts = 5 (dashed line) and 
rs — 10 (dot-dashed line) caused by interactions, is the 
radius, in units of the Bohr radius, of a sphere containing 
one electron and in metallic elements at atmospheric pressure 
ranges from roughly 2 to 6. The unit discontinuity of the 
momentum distribution for the noninteracting system (solid 
line) at the Fermi surface is reduced by 28% and 42% respec- 
tively. Corresponding increases in the kinetic energy over the 
noninteracting case are 30% and 60%. Bottom figure: For 
''He at T = and P ^ roughly 92% of the zero momentum 
condensate (black vertical bar along y axis) in the noninter- 
acting Bose system is promoted to higher momentum states 
(distribution given by dashed line) leading to a kinetic energy 
increase of over 14 K per particle. 

This argument fails however, when generalized to 
nonzero temperature, T, where the system occupies a 
thermal distribution of energy eigenstates. From the 
Gibbs variational principle, the free energy functional, 

F[p] - TT[Hap] + keT Tr[p In p] , (2) 

takes its minimum value F[po] — Kq — TSq for the equi- 
librium density operator po = e~''^''/Tr[e~'^^''] where 



2 



f3 = l/ksT and fc^ is Boltzmann's constant. Any other 
normalized density operator, such as that for the inter- 
acting system, p — e^^^ /Ti[e~^^], gives a higher value, 

F[p]^K~TS > F[po]^Ko~TSo. (3) 

The ground state inequality, Eq. |l|, now generalizes to, 

K > Ko + T{S - So) . (4) 

Since interactions often increase order, the entropy S 
of the interacting system can be less than that of the 
noninteracting, Sq, and S — Sq may be negative. This 
allows Eq. ^ to be satisfied while the ground state in- 
equality, K > Kq, is not. In general, at nonzero tem- 
perature nothing forbids interactions from lowering the 
kinetic energy. 

This kinetic energy lowering would only be observed 
in an intermediate range of temperatures since in the 
high temperature, classical limit, interactions only effect 
equilibration rates, not the final Maxwellian momentum 
distribution. 
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FIG. 2: Behavior of free fermion (top figure) and free boson 
(lower figure) radial distribution functions for indicated tem- 
peratures measured in terms of tlie Fermi temperature Tf or 
condensation temperature Tc respectively. 



WEAK COUPLING LIMIT 

Not only is this lowering theoretically possible, it can 
be shown to occur for a variety of weakly interacting 
physical systems using lowest order thermodynamic per- 
turbation theory. The change in the Helmholtz free en- 
ergy, F, to first order in the interaction potential is 



F^Fn 



{U)o 



(5) 



where {U)q is the potential averaged over the free particle 
configurations. Since the total energy is given by, E — 
d{PF)/dp, the first order change in the kinetic energy 
from Eq. ^ is 



dp 



(6) 



For a single component system of N particles with pair 
interaction, V{r), the first order energy change reads. 



(7) 



where n is number density. The first order kinetic energy 
change 
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dT 



V{r) dr^ 



(8) 



depends on the temperature dependence of the free par- 
ticle radial distribution function, go(r). 

As shown in Fig. |^, for both fermions and bosons 
dgQ{r)/dT can be positive at all r leading to a kinetic 
energy decrease for repulsive interactions although for 



quite different reasons in the two cases. For fermions the 
increase of go{r) with temperature is due to the filling 
in of the "Fermi hole" , the region near the origin where 
due to the Pauli principle like spin particles are excluded. 
The size of this region is roughly the de Broglie thermal 
wavelength which is proportional to I/VT. 

For noninteracting bosons at T = 0, all particles are 
in the zero momentum condensate and go{r) — 1. As 
the temperature increases particles are promoted from 
the condensate and ga{r) increases at all r for T/Tc < 
(2/5)^/"^, where Tc is the Bose condensation temperature, 
and thereafter at small r, reaching go{0) = 2 at T = Tc ||, 
H . Since Tdgo (r ) / dT vanishes at T = , the ground state 
inequality, Eq. |l|, is not violated. 

The mechanism for the narrowing of the electron gas 
momentum distribution can be seen from the first order 
shift in the free electron energy levels M, 



Ae(k) = - 
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"o(k) 

k-k'p 



(9) 



The decrease of no(k) with k leads to a larger lowering 
of the energy levels at smaller k and thus to a population 
redistribution and momentum distribution narrowing as 
demonstrated in Fig. || 

SIMULATION RESULTS 

For the electron gas, this weak coupling expansion ap- 
plies only at very high densities, - 0.5, and is of httle 
use for realistic condensed matter densities. It is neces- 
sary to use more powerful quantum many body simula- 
tion methods. Two such methods, path integral Monte 
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FIG. 3: Example of momentum distribution narrowing in 
the fully polarized electron gas at Ts = 4 and T — Tpormi from 
simulations with 57 particles. 



FIG. 4: Relative excess kinetic energy vs temperature for the 
unpolarized electron gas at indicated values as calculated 
from PIMC simulations with 66 particles in the periodic cell. 



Carlo (PIMC) § | and a variational trial density ma- 
trix method (VDM) have been used here to search 
for narrowing of the momentum distribution in a variety 
of fermion systems. 

Liquid '^He is a strongly coupled system. Due to its 
steeply repulsive short ranged interactions the weak cou- 
pling arguments of the preceding paragraphs do not ap- 
ply. Restricted PIMC calculations using free particle 
nodes for '^He at number density n — 0.01636A and 
temperatures from 1 to 40 K show the kinetic energy to 
be more than 5 K above the noninteracting value. Nar- 
rowing of the momentum distribution is not found for 
this system. 

For the electron gas model both restricted PIMC and 
VDM calculations find momentum distribution narrow- 
ing. The effect is small, with the relative decrease of the 
kinetic energy from its ideal value less than a few per- 
cent |lO[ |ll[, but well within the accuracy of the two 
methods. The computation of momentum distribution 
for fermions with PIMC (see Fig. ||) required extending 
the restricted path integral technique to simulations with 
open paths jl^, |l^ . 

Fig. ^ shows the difference in kinetic energy of the 
homogeneous electron gas for different densities and 
temperatures indicated by the degeneracy parameter 
T/Tpermi- For densities corresponding to rs ~ 4, a low- 
ering of the kinetic energy with respect to the noninter- 
acting value is found. The magnitude of the lowering 
increases with density (smaller r^). The statistical un- 
certainty of these results was estimated at = 1 and 
T — Tpcrmi by a study of finite size effects (using 14, 
38, 66, and 114 particles in the periodic simulation cell) 
and path discretization errors (using 16, 32, and 64 time 



steps) which gave {K - Ko)/Ko = 0.023 ± 0.002. The 
region of the effect, as predicted by PIMC simulations, 
is shown in the high temperature and density phase dia- 
gram in Fig. H 

Simulations of dense hydrogen plasma have also 
indicated kinetic energy lowering, which we confirmed 
and extended. We found a maximal lowering of (K — 
Kq)/Ko = 0.007 ± 0.001 for r,. = 0.5 and T = Tpcrmi- 
The magnitude is reduced compared to the electron gas 
because Coulomb interactions of electrons and protons 
counteract the entropic lowering. This also results in 
a reduction of the parameter region in the temperature 
and density plane for which the effect can be observed. 
However, in the limit of high density, the boundaries 
for the electron gas model and for hydrogen converge 
because interactions become weaker in the high density 
limit. Fig. |^ shows that the lowering region of hydrogen 
includes conditions near the core of the sun or other low 
mass stars as well as on the compression path of in- 
ertial confinement fusion ||l5|, indicating that hydrogen 
is a potential candidate for experimental verification of 
this effect. 



CONCLUSIONS 

Analytical and numerical arguments show that interac- 
tions can reduce the kinetic energy of a quantum system 
below the corresponding noninteracting value. The asso- 
ciated narrowing in the momentum distribution can po- 
tentially be verified experimentally. Prime examples are 
the Bose condensates in magnetic traps, in which the mo- 
mentum distribution can be inferred by measuring the ex- 
pansion rates after the magnetic field has been switched 
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FIG. 5: Proposed region for momentum narrowing {K < Ko) 
in the unpolarized electron gas (area above the sohd line) 
and in hydrogen plasma (area above the dot dashed line) as 
derived from PIMC simulations. The dashed line indicates 
the compression path of inertial confinement fusion |^| . The 
conditions near the core of the sun and a low mass star (0.3 
solar masses) are also indicated 

off. Sufficiently accurate measurements of the momen- 
tum distribution could be used to provide information 
about the nature of the interactions of the trapped atoms. 

Lowering of the kinetic energy was demonstrated in the 
homogeneous electron gas model, one of the fundamental 
models in condensed matter physics and a well studied 
example of a one-component plasma. Corrections to the 
noninteracting kinetic energy are always positive at zero 
temperature and vanish in the high temperature limit, 
which makes it counterintuitive to expect lowering at in- 
termediate temperature. Research on the one-component 
plasma has focused more on the internal energy, pressure 
and correlation energy. However, the kinetic energy is 
relevant of our understanding of quantum systems and 
can be determined experimentally by measuring the mo- 
mentum distribution. 

The effect, previously indicated but left uninterpreted 
in dense hydrogen p^ , is reconfirmed here. The argu- 
ments presented above suggest that the kinetic energy 
reduction should be present in weakly coupled systems 
with repulsive pair interactions. Simulation evidence for 
this effect was presented here only for Coulombic sys- 



tems. More work is necessary to understand how these 
arguments apply to other systems and where the entropic 
reordering in the thermal population of states leading to 
a lower kinetic energy can best be observed experimen- 
tally 
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